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Abstract 

A heat kernel approach is proposed for the development of a general, flexible, and 
mathematically tractable asset pricing framework in finite time. The pricing kernel, 

r— i giving rise to the price system in an incomplete market, is modelled by weighted heat 

kernels that are driven by multivariate Markov processes and that provide enough de- 

PLh grees of freedom in order to calibrate to relevant data, e.g. to the term structure of 

bond prices. It is shown how, for a class of models, the prices of bonds, caplets, and 
swaptions can be computed in closed form. The dynamical equations for the price 
processes are derived, and explicit formulae are obtained for the short rate of interest, 
the risk premium, and for the stochastic volatility of prices. Several of the closed-form 
asset price models presented in this paper are driven by combinations of Markovian 

t-h jump processes with different probability laws. Such models provide a rich basis for 

consistent applications in several sectors of a financial market including equity, fixed- 
income, commodities, and insurance. The flexible, multidimensional and multivariate 

qq structure, on which the asset price models are constructed, lends itself well to the trans- 

parent modelling of dependence across asset classes. As an illustration, the impact on 
prices by spiralling debt, a typical feature of a financial crisis, is modelled explicitly, 
and contagion effects are readily observed in the dynamics of asset returns. 

Keywords: Asset pricing, pricing kernel, Markov processes, Levy random bridges, 
weighted heat kernels, equity, interest rate derivatives, debt, spread dynamics, and 
^ contagion. 

Introduction 2 

Pricing kernel models and the pricing of bonds, caplets, and swaptions 4 

Closed-form and explicit price models 7 

Dynamical equations 10 

Incomplete market models driven by LRBs 12 

General asset pricing in finite time 18 

Spiralling debt and its impact on global bond markets 22 

References 26 



1 



1 Introduction 



In this paper, we shall take the view that in a modern asset pricing framework (i) pricing 
models should be coherent across all asset classes traded in a financial market, (ii) securities 
pricing used in the front offices of financial firms should be compatible with asset risk 
management, and (iii) pricing formulae should be applicable in the banking industry and 
also in the insurance sector. Expressed in other words, these three requirements state that 
modern pricing models ought to be consistent under the real probability measure P and 
the risk-neutral measure Q, and at the same time they should retain a high degree of 
flexibility and mathematical ease while guaranteeing the coherence of the price system for 
all financial assets. 

In what follows, we propose an asset pricing framework that can be applied, in princi- 
ple, to all asset classes and that is mathematically tractable so that Monte Carlo techniques 
are never necessary for scenario simulations of asset price dynamics. The proposed ap- 
proach includes partial automatic calibration to market data such as initial prices of assets. 
The price system of assets traded in a financial market shall be developed by modelling the 
pricing kernel (state-price density) first. Once the stochastic framework for the pricing ker- 
nel is built and the connection with bond prices is established, we go on to show how price 
processes for other asset classes can be derived in a natural way. We also consider how the 
situation, in which the debt of a sovereign country gets out of control, can be incorporated 
in the same pricing framework without introducing extra assumptions to include effects of 
credit risk. 

The general setup of the asset pricing framework is developed in finite time, t e [0, U] 
for U < oo. We model a financial market by a filtered probability space (f2, J^,P, {J^J) 
where P denotes the real probability measure, and where {& t } is the market filtration. 
We consider a (multi-dimensional) process {X t } on (f2, J^,P), and assume that the market 
filtration is generated by {X t }. Furthermore, we assume that {X t } has the Markov property 
with respect to {& t }> its natural filtration. Then, we introduce the pricing kernel process 
{7i J to model the market agent's preferences and the dynamics of interest rates in the 
economy which {n t } is associated with. We write {S t } < t < T<u for the price process of a 
dividend-paying asset, and let {D t } < t < T<u denote the (continuous) dividend stream up 
until T. Then the price S t at time t is given by 



In order to calculate asset prices explicitly, the following ingredients need to be specified: 
(i) The Markov process {X t } that generates the market filtration, and thus the market in- 
formation; (ii) the pricing kernel {n t }, and thus the dynamics of the interest rates and 
the agent's preferences; (iii) the random variable S T and the process {D t }, thus the asset's 
terminal cash flow and the dividend stream, respectively. All ingredients are specified in 
such a way that the price process {S t } is adapted to the market filtration generated by {X t }. 
For textbooks about the theory of pricing kernels, preferences, asset pricing, and interest 
rates modelling, one may consult, e.g., Back (2010), Bjork (2009), Cochrane (2005), Duffie 
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(2001), and Brigo & Mercurio (2006). 

In the next section, we introduce weighted heat kernels to define the class of pricing 
kernels, and thus give rise to the asset pricing framework proposed in this paper. Heat ker- 
nel models for the development of stochastic price systems have been proposed by Akahori 
et al. (2009) in an infinite-time setting, and more recently by Akahori & Macrina (2012) 
in a finite-time context. We further develop the weighted heat kernel pricing models intro- 
duced in the latter work, and derive formulae for the price processes of discount bonds and 
swaptions, and the associated interest rate process. The stochastic short rate of interest is 
by construction non-negative. 

In Section [3] we construct a class of pricing kernels which leads to closed-form price 
processes for bonds, caplets, and swaptions — all in one. Explicit price models are then 
obtained by specifying the dynamics of the market information flow and the degrees of 
freedom in the formulation of a particular pricing kernel model. The dynamics of certain 
price processes have time-dependent lower and upper bounds, a feature we not necessarily 
view as being a shortcoming. 

In Section |4| we derive the dynamical equation for the bond price processes introduced 
in the previous section for the case that the market filtration is generated by a diffusion 
process. The market price of risk process is also obtained endogenously which, for this class 
of pricing models, incorporates a discernible part that can be identified as the incentive for 
accepting model risk. 

In Section [5] we introduce multivariate Levy random bridges and extend the pricing 
framework to an incomplete market. A formula for the computation of option prices by 
Fast-Fourier-Transform methods is given, and asset price models with higher-order rational 
structures are developed. 

In Section [6] the proposed pricing kernel approach is applied to general asset pricing, 
and we show how asset price models constructed under the P-measure translate into as- 
set price models equipped with stochastic interest rate and stochastic volatility under the 
Q-measure. The interaction between the bond price process and, e.g., the equity compo- 
nent of the discounted share price process is clearly identifiable. This property renders the 
herewith proposed asset pricing framework also appealing for the construction of hedging 
strategies against losses due to the exposure of a financial position to a specific market 
sector. Additional examples of explicit multivariate asset pricing models driven by jump 
processes are given. 

In the last chapter, spiralling sovereign deficit is modelled and its impact on the price 
dynamics of sovereign bonds is shown. The flexibility of the considered heat kernel state- 
price density approach allows for the construction of explicit dependence models linking 
the price evolution of, e.g., bonds issued by several sovereign governments. Contagion 
effects arise endogenously, and the graphs in Section [7] illustrate the impact of dependent 
economies and markets on the price dynamics of assets. 

The list of investigations connected with the asset pricing approach presented in this pa- 
per is by no means complete. The developed theory is applicable to the modelling of foreign 
exchange rates and the pricing of foreign exchange securities. Inflation-linked bonds and 
other indexed assets can be priced in a similar way. Then, there are commodity assets, in- 
cluding energy and agricultural products, which may involve insurance contracts to hedge 
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against substantial losses due to averse weather conditions. In the last section of this paper, 
effects on asset prices by perceived solvency risk are considered in the dynamics of their 
returns. Further work in this direction, may include the modelling of "costs of funding" 
arising from shifts in asset prices due to the deterioration of a creditor's economic situation. 
Even though full-fledged credit risk models have not yet been developed within the present 
asset pricing framework (and this is yet another project), the material in Section[7]may also 
be looked at from the perspective of modelling "credit valuation adjustments". The bounds, 
within which the dynamics of certain price processes are confined, could be exploited to 
model the levels of sustainability for the costs of funding, or applied by financial regula- 
tors to impose time-dependent capital requirements. Although treated in this paper, the 
explicit modelling of dependence structures for asset portfolios remains somewhat in the 
background. The weighted heat kernel approach however offers a versatile basis for the 
inclusion of manageable dependence models that could be useful for the risk analysis of 
portfolios. One could begin with mean-variance optimization whereby covariance matrices 
are explicitly modelled, and partial information about the underlying risk factors drives the 
optimal portfolio. Another investigation may concern models for volatility surfaces, which 
arise from the selection of particular pricing kernels and their application to specific asset 
classes. Such an investigation extends to the analysis of the derived option price models 
and their calibration to data relevant for the pricing of assets. 



2 Pricing kernel models and the pricing of bonds, caplets 
and swaptions 

We proceed to specify the form of the pricing kernel model {n t }. We consider the following 
class: 



r u-t 



rc t =/o(0+/i(0 



E[F(t + u,X t+u )\X t ]w(t,u)du, (2.1) 



where {X t } is an unspecified Markov process defined for t e [0, U] such that t+u < U, / (t) 
and fi(t) are deterministic, non-negative, non-increasing functions, F(t,x) is a positive 
measurable function, and w(t,u) is a so-called weight function satisfying 

w(t,u -5) < w(t -s,u) (2.2) 

for s < t A u. We introduce the following short-hand notation: 

p(t, u,X t ) := E [F(t + u,X t+u ) \X t ] . (2.3) 

It can be proven that the considered pricing kernel processes are indeed supermartingales 
adapted to the filtration generated by {X t }. We refer to Akahori & Macrina (2012) for a 
proof that can be applied also in the present context. 

As a special case of the general pricing formula ( |1.1| ), the price process of a discount 
bond with maturity T, is given by 

1 r 

P tT = —E[n T \X t ], (2.4) 
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where < r < T < U. We keep in mind that the market filtration is generated by the 
Markov process {X t }, and thus it suffices to take the expectation conditional only onX t . 
The conditional expectation of ti t can be computed explicitly to obtain 



r u-t 



E[n T \x t ] =/ (r)+/ 1 (r) 



p{t,u,X t )w{T,u- T + t)du, 



(2.5) 



T-t 



where the tower property is invoked and a variable substitution is applied. We define 



r u-t 



Y tT = 



p{t,u,X t )w{T,u- T + t)du. 



(2.6) 



T-t 



The bond price process can then be written in the compact form 

fo(T) + A(T)Y tT 



P tT = 



/o(0 + /i(t)r tt ' 



and the initial term structure is given by 



Pot 



/o(0+/i(Oro t 

/o(0) + /i(0)Y c 



oo 



We deduce that 



(2.7) 



(2.8) 



(2.9) 



/o(0 = P 0t [i+/i(o)y 00 ] -fi(t)Y ot , 

where we may set f (0) = 1 with no loss of generality. By inserting ( 2.9 ) in ( 2.7 1, we obtain 

PoT+YiT) (Y tT -Y 0T ) 



where 



?tT Pot+y(t)(Y tt -Y ot ) ' 

m /1(t) 



(2.10) 



for < t < T. Similarly, the expression for the pricing kernel can be written in the form 



(2.11) 



K t = n [p ot +y(t) (y„-y ot )], 



(2.12) 



where 7i = 1 +/i(0)Yo . Assuming that the bond price function is differentiable with 
respect to T, the expression for the instantaneous forward rate {r tT } is given by 



"tT = ~ d T ln {Pit) , 

d T P 0T + (Y tT - Y 0T ) d T y{T) + y(T) (d T Y t 



P 0T + y(T)(Y tT -Y 0T ) 
The process {r t } for the short rate of interest can then be deduced by setting r t = r tT 



(2.13) 



T=t 



r t = 



djPot + {Y tt - Y 0t ) {d T y{T)} T =t + y(0 { ( d T Y tT ~ d T Y QT ) } T=t 
Pot + y(t)(Y tt -Y 0t ) 



(2.14) 
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Once the bond price system is derived, one can calculate the price of fixed-income 
derivatives such as caplets and swaptions. We consider a t-maturity swaption contract with 
strike K, which is written on a collection of discount bonds P tT . with maturities {T^i n . 
An application of the pricing formula shows that the price Swp Qt of the swaption at 
time zero is 

nJl-P^-K^P^) . (2.15) 



Swp 0t = — E 

7l 



i=l 



where 



P t T = 



P 0Ti +y(T i )(Y tTi -Y 0Ti ) 
Pot + y(t)(Y tt -Y 0t ) 



Then, by use of (|2.12|), we obtain 

P ot -P 0Tn -Kj]p 0Ti +y(t) (Y tt - Y ot ) 



Swp 0t = E 



(2.16) 



(2.17) 



-y(T n ) (Y tTn - Y 0Tn ) -K^yiTj (Y tT{ - Y 0T .) 



i=l J 

The price of caplets can be calculated in an analogous way. Further details for the calcula- 
tion of caplets and swaptions follow in the next section. 

We may wonder at this stage whether it might be possible to construct a class of discount 
bond price processes for which the associated prices of interest rate derivatives can be 
calculated in closed form. We shall present such bond price models in detail in the next 
section. Before, we prepare the ground by making the following observation. 

Proposition 2.1. Let {M t } < t<(7 be an {& t }-adapted F-martingale that induces a change-of- 
measure from P to an equivalent auxiliary probability measure M. Analogous to (2.6 ), let 

i Y ^}o<t<T<u be defined by 

r u-t 



v m _ 



p M (t, u,X t ) w(T,u-T + t) ciu, 



T-t 



where 

Then the process 



p M (t,u,X t ) = E M [F(t + u,X t+u )\X t ]. 



n t = n [P 0t + y(0 (Y™-Y™)]M t 
is a positive {{& t } ,F)-supermartingale where 

/o(t)+A(t)y* 



Pot 



and 



y(0 = 



/o(o)+/i(o)y*f : 
hit) 



l+/i(0)F. 



(2.18) 

(2.19) 
(2.20) 

(2.21) 

(2.22) 



oo 
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Proof. By construction {n t /M t } is a positive ({J 2 ^},] 
< t < u < U, we have 



-supermartingale. That is, for 



^[n u \X t ] = ^[n [p 0u + y{u)(Y™-Y™)]M u \X t 
= ^[n [p 0u+ y(u)(YZ-Y^]x t ]M t , 
< n [P 0t +y(t)(Y*-Y*)]M t = 7i t . 



(2.23) 



□ 



Thus, if convenient, we can construct pricing kernel models by considering M-propagators 
{p M [t,u,X t )} while using the form (2.20). In such a case, the form of the asset price 
formulae remain unchanged, and for the discount bond price process one has 



PtT = 



P T+y(T)(Y™-Y™) 

p ot +y(t)(Y™-Y™) ' 



(2.24) 



3 Closed-form and explicit price models 

In this section, we construct a class of pricing kernels for which the price processes of 
underlying and derivative assets are obtained analytically. We explicitly calculate the price 
processes of bonds, caplets, and swaptions, and note that the derived examples lead to 
bounded price processes — a property we not necessarily view as a shortcoming. 

Proposition 3.1. Let {M t } < t<u define a measure M on & t such that M is equivalent to F 
onte [0, [/). Let {A t } < t<u be an {{& t } ,Wf)-martingale and let b(t) be a deterministic and 
non-increasing function such that n t > 0/or all t e [0, U). Then, 



n t = [-Pot + Ht)A t ] M t 
is a pricing kernel, and the discount bond system P tT takes the form 

P 0T + fa(T)A £ 



P tT = 



P 0t + b(t)A t 



(3.1) 



(3.2) 



Proof. To verify that (3.1 1 is indeed a positive supermartingale, one follows the proof 



of Proposition 2.1 The discount bond price process (3.2) is derived by applying (|2.4[) and 



by observing that 



Pa 



E F [[P 0T + b(T)A T ] M T \X t ] 

[Pot + Kt)A t ] M t 
P 0T + b(T)E M [A T \X t ] 



P 0t + Kt)A t 

Here we have used the Bayes formula to change the measure from P to 



definition an ({^ t },M)-martingale, the result (3.2) holds. 



(3.3) 

Since {A t } is by 
□ 
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Lemma 3.1. The short rate of interest process {r t } associated with the pricing kernel models 

d t p ot +A t d t b(t) 



(3.1) is of the form 



r f = -- 



P 0t + b(t)A t 



(3.4) 



Proof. We first derive the instantaneous forward rate {r tT } as in (2.13), and then we 

□ 

r=t 



set r t = r 



tr 



Lemma 3.2. In the pricing system specified in Proposition (3.1), the prices Cp 0t and Swp 0t 
at time zero of t-maturity caplet and swaption contracts are as follows: 



Caplet: 

C Pot = -{d t P 0t + KP 0t ) 
where P[A t e da] = p(a)da and 



ap(a)da- [d t b(t) + Kb(t)] 



ap(a)da, 



(3.5) 



a : a < — - 



d t P 0t +KP t 



Ot 



d t b(t) + Kb(t) 



(3.6) 



Swaption: 
Swp 0t 



Pot -Pi 



or,, 



■K 



1=1 



OT, 



p(a) da 



+ 



i=l 



ap(a)da, (3.7) 



where P[A t e da] = p(a)da and 



a := < a : a > 



^2;=i ^or, — Pot +-^or„ 



(3.8) 



b(t)-b(r n )-K^ =1 b(rO j' 

Explicit models. In order to obtain explicit pricing models, the following quantities 
need to be specified in the definition of the pricing kernel (2.1 1: (i) The finite-time Markov 
process {X t } that generates the market filtration and drives all prices, (ii) the positive func- 
tion F(t + u, x) that, to a great extent, characterises the type of pricing model, (iii) the 
weight function w(t,u), and (iv) the deterministic functions / (t) and/ 1 (t). The particular 
class of models considered in this paper allows for explicit calibration of / (t) to the initial 
term structure, and for a one-to-one correspondence between the degree of freedom /j(t) 
and option data (e.g., caplets and swaptions). Thus we only specify the Markov process 
{X t }, F(t + u, x), and w{t,u). In the following two examples, we let the Markov process 
{X t } be given by the Brownian random bridge {L tU } defined by 



L tU — a tX v + [3 tU , 



(3.9) 
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where a is a constant parameter, X v is a random variable with a priori density p(x), and 
{Ptu}o<t<u i s an independent standard Brownian bridge where U is fixed. The change- 
of-measure density martingale {M t } in (3.1) that is associated with such an information 
model, satisfies 

dM t = - ^^E p [Xu | L tu ] M t dW t F (3.10) 



where {W t } is defined by (4.6) in the next section. 

Quadratic models: We choose F(t + u, x) = x 2 , and set w(t,u) = U — t — u. For this class 
of models, we obtain 



U 



(U-t) 2 tu U-t 



and 



b(t) = 



(c/-0 4 /i(0 

4[/[l + ^/ 1 (0)[/ 3 ]' 



(3.11) 



which determine explicitly the processes for the pricing kernel, the bond price, and the 
associated interest rate. Furthermore, we can now work out the explicit expressions for the 
caplet and swaption prices via ( |3.5| ) and ( |3.7| ), respectively. For the caplet price, we obtain 

1 {d t b(t)+Kb(t))-(d t P 0t +KP 0t ) 



Cp 



0t 



U-t 



N(k) 



+ [8 t b(t) + Kb(t)] -^expf-l-K 2 ), (3.12) 

U — t v27i V A J 



where iV(x) denotes the cumulative normal distribution function, and k is defined by 

(3.13) 



U-t 



d t p ot +KP 0t 



U-t d t b{t) + Kb{t)^ 

Exponential quadratic models: In this case we choose a function F(t + u,x) that depends 
explicitly on time, that is: 

F(t+ "' x) = exp ( 2 a/- t - U ))' 



(3.14) 



and take the weight function to be w(t,u) — (U — t — u) v l ^ 2 where 17 > 1/2. Then, we 
have: 



A t = 



U 



exp 



J tu 



2{U-t) 



1 and b(t) = 



{U -t^-^u 112 /^) 

l+/i(0)[/" ' 



(3.15) 



As for the quadratic models, the pricing kernel, the bond price, and the associated interest 
rate are all determined in closed form. For the price of a swaption contract, we obtain, for 
this class of models, the following: 



Swp 0t 



P 0t -P 0Tn -Kj]P 0Ti )N(-v) 



i=l 



+ 



b(t)-b(T„)-ic£|b(rO 



i=l 



N(v) -N (v \fl - t/L/j 



, (3.16) 
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where 



v := 



\ 



2U 



In 



- 1 / 2 /" ^Sr=i p o ri + p OT n -P. 



Ot 



+ 1 



(3.17) 



Analogous applications of Lemma 3.2 lead to the explicit expressions for the price of a 



swaption utilising the quadratic models, and for the price of caplets in the case of the ex- 
ponential quaratic models. Although written in a less unifying form — and in the case of the 
quadratic models, including less useful degrees of freedom — the quadratic and exponential 
quadratic models were first developed in Akahori & Macrina (2012). We also note that 
interest rate models with a quadratic or exponential quadratic structure have been studied 
in Jamashidian (1996) and McCloud (2009, 2012), too. 



Boundedness of prices. Bond prices fluctuate by construction between zero and one, 
and the associated interest rate is non-negative. However, the bond price processes pro- 
duced by the above models have tighter bounds, and the same holds for the interest rate 
and the yield of the bond. As a consequence, the strike price of call bond options have 
an upper bound for in-the-money options. One might think that having bounded bond 
prices and associated interest rates is a shortcoming. On the contrary, we think that such 
a feature may be advantageous, especially if the time-dependent bounds are wide enough 
for the interest rate to have sufficient freedom. The bounds may be put in relation with 
economic policies of which goal is to keep bond prices within a certain range. In turn, this 
may suggest to use the additional degree of freedom /i(t), cast inside the deterministic 
function b(t), to include regulators' policies, for instance. Research regarding bounded as- 
set prices and the relation to regulators' policies and markets shall be continued elsewhere. 
We shall keep the boundedness property inherent in certain rational asset pricing models 
in our mind for when we later turn to general asset pricing, Section [6j and to the impact on 
prices by an economy's spiralling deficit, Section [7j 



4 Dynamical equations 

We derive the stochastic differential equation of the bond price for the case that the mar- 



tingale {A t } in Proposition 3.1 is a continuous diffusion process. The drift term of the 
stochastic differential equation reveals the market price of risk, which is obtained endoge- 
nously. It turns out that the market price of risk process is constituted by two distinct 
components. One part is associated with the stochasticity of a financial market due to noisy 
information about the market factors. The second component of the market price of risk 
can be identified as model risk that is directly related to the choice of the class of price 
models. Furthermore, the Brownian motion that drives the bond price process emerges 
also endogenously and in fact it is the innovations process updating the price process as 
the quality of market information improves. 
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Proposition 4.1. Let {W f F } be a standard ({& t },F)-Brownian motion. Let the (|^ t },M)- 
martingale {A t } satisfy 

dA t = v t (dW t p + # t dt), (4.1) 

where {v t } and {# t } are well-defined {& t }-adapted processes. Then the dynamical equation of 
the discount bond with price process (3.2) is given by 

dP tT 



tT 



where 



(r t + ?i t n tT )dt + n tT ciW^, 



d t p Qt +A t d t b{t) 



n tT = ■ 



Pot + b(t)A t 



P 0t + b(t)A t ' 
6(T) 6(t) 



P 0T + b(T)A t P ot + 6(t)A £ 



Proof. We apply Ito's Lemma to (3.2) and use (4.1 ). 



(4.2) 

(4.3) 
(4.4) 
(4.5) 
□ 



By specifying the driving Markov process and by selecting a particular type of process 
{A t } in Proposition (4.1 ), we also obtain the result in the following lemma. 

Lemma 4.1 . Let the Markov process generating the market filtration {& t } be the information 
process (3.9). Then the Brownian motion {W f F } satisfies 



dW, 



t = dl tp - — {crUE v [X v I L tu ] - L tu ) dt, 



and the process {6 t } in Proposition 4. 1 is given by 

oU 



U-t 



E F [X v | L tu ] . 



The process {v t } in Proposition 4.1 satisfies 



2U 



(u-ty 



in the case of the quadratic models (3.11), and it satisfies 



u-t A 1/2 



v f = 



u-t 



Ltu exp 



J tu 



2(17 -t) 



(4.6) 



(4.7) 



(4.8) 



(4.9) 



in the case of the exponential quadratic models ( 3.15). 
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The process {W t M }, satisfying dW f M = dW t P + # t dt, is an ({J r J,M)-Brownian motion. 
This follows from the Girsanov theorem. We emphasize here that the risk premium process 
{X t } is constituted by the process {tfj and {v t }. We observe (a) that {~& t } is determined 
by the filtration model, that is, by the choice of the generating Markov process {L tU }, and 
(b) that {v t } depends on the selection of the heat kernel models, that is, on the choice for 
F(t, x) and w(t,u). We can view {# t } as the risk premium component associated with the 
uncertainty in the market modelled via the information flow process {L tU }. The component 
{v c } may however be interpreted as the premium associated with model risk since it is 
closely related to the choice of the specific asset price model. By applying Ito's formula to 



(2.13 ), we obtain the following result: 



Proposition 4.2. The instantaneous forward rate {r tT } of the bond price process (3.2) satisfies 
the dynamical equation 

dr 

— — = -cr tT n tT dt + cr tT (dwf + A t dr) , (4.10) 

r tT 

where the instantaneous forward rate volatility {<J tT } is defined by 

~ Vt [pot + KT)A t ~ d T P 0T +A t d T b(T)J • (4 - n) 



Applying Ito's formula to (2.141 gives the stochastic differential equation for the short 



rate of interest {r t } in the same set-up: 



dr 

— = /x t dt + o- t (dW t P + A t dt), (4.12) 



where {A t } is the instantaneous market price of risk (4.41, cr t = 0"tT\ T _ > an< ^ 



= d t p 0t +A t d t b(t) _ d n p 0t +A t d tt b(t) 
Mt p 0t + b(t)A t d t p ot +A t d t b(t) ■ ■ 

The processes driving the dynamics of the presented models have Gaussian laws, and 
the dynamics of the deduced instantaneous forward rates have the HJM-form, c.f. Heath 
et al. (1992), Filipovic (2009). It can be shown that {W t F } < t<!7 is an ({J r J,P)-Brownian 
motion, see, e.g., Brody et al. (2008). 



5 Incomplete market models driven by LRBs 

In this section, we extend the pricing framework to include multi-dimensional risk factors, 
and we generate asset pricing models in an incomplete market. We consider a class of 
finite-time Markov processes, the so-called "Levy random bridges" (LRBs), as constructed 
in Hoyle et al. (2011). An LRB can be interpreted as a Levy process that is bound to 
have a prescribed, albeit arbitrary, distribution at a fixed future time. The LRB and the 
generating Levy process are linked by a an equivalent probability measure with respect to 
which the LRB has the law of the generating Levy process. Before we slightly extend this 
result appearing in Hoyle et al. (2011), we first give the definition of a multivariate LRB: 
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Definition 5.1. We say that {L t [/}o<t<[/ is a multivariate LRB on R m if the following are 
satisfied: 

1. The random variable L uv on R m has multivariate marginal law v{z\ z e R m . 

2. There exist a multivariate Levy process {L t } <t<[/ on ^ m sucn that L t has multivariate 
density function p t (x) on R m for all t e (0, 17]. 

3. The marginal law v(z) concentrates mass where p v {z) is positive and finite, that is 
< Pu(z) < oo for v(z) almost every z. 

4. For every n e N + , every < t 1 < . . . < t n < U, every (x 1; . . . , x n ) e R m x R" and v(z) 
almost every z, we have 

¥[L tiU <x 1 ,...,L tnU <x n \L uu = z] =¥[L ti <x 1 ,...L tn <x n \L u = z]. (5.1) 

Proposition 5.1. Let {L t u}o<t<u denote a multivariate LRB with marginal law F[Luu e 
dz] = v(dz). Let the multivariate Levy process {L t } < t < v , which generates the LRB, have 
density p t (x)/or all t e (0, 17]. t/nder the measure L defined by 



' c " dL 



— v(dz), (5.2) 



the L7?B {L tU } has the law of the generating Levy process for t e [0, U). 

Proof. The verification of this proposition follows closely the results leading to Proposi- 
tion 3.7 in Hoyle et al. (2011). □ 

The measure L is rather useful for several calculations. LRBs, which have joint marginal 
law at t = U and which are generated by independent Levy processes, are independent 
under L: 

Proposition 5.2. Let {L^o^y be a multivariate LRB on R m generated by a Levy process 
on M. m of which multivariate density function p t (z) factorises, that is 

m 

p t (z) = Y\p*(z k ). (5.3) 

fc=i 

Under L, the components, {L^} and {L^} j, of the LRB on R m are independent and each 
LRB component has the law of the respective component of the generating Levy process {L t }, 
for t e [0, [/). 

Proof. We show that the generating function of the multivariate LRB factorises under 
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L. We have: 



E 



= E L 



exp 



.(fc) 
tu 



\k=l J 
exp M>] a fc y k J L 

exp J]a fc y fc ]~[L [l^ 3 S dy fc 



L£edy 1 ,...,L[* ) edy n 



.fc=i 



fc=i 



m f 

= n exp y ^ l [ L tu e ^ 

fc=i 

m 



fc=i 



(5.4) 



□ 



Next we make use of the L-independence property of the LRBs to propose multi-factor 
pricing kernel models in the situation where the driving Markov process is a multivariate 
LRB. We assume that the market filtration {& t } is generated by a two-dimensional LRB, of 
which first component is a Brownian random bridge, 



■CD 



aX™ t + p tu , 



(5.5) 



and of which second component is a gamma random bridge defined by 

4?=4 2) r tt /- (5.6) 

The Brownian bridge {I3 tu } and the gamma bridge {y tU } are assumed independent of each 
other and also independent of the random variables and X^\ However, the two X 
random variables are dependent and have a priori joint marginal law v(z 1 ,z 2 ). This set-up 
is of the kind considered in Proposition 5.2 Next we recall the bond pricing formula (2.24 1 
where, this time, 



r L = 



r u-t 



T-t 



r U-t 



w{T,u x -T + t,u 2 -T + t) 



T-t 



x E 



F[t + u 1 ,t + u 2 , L^ Ui>u , lfl U2>u ) | ig, L%] & Ul du 2 . (5.7) 



Since, for t e [0, U), the two LRB components have each the law of the corresponding 
underlying Levy process component — that is respectively the Brownian motion and the 
gamma process — the conditional expectation simplifies considerably under L. That is, 



E- 



F (t + u a , t + u 2 , L^ UuU , lfl 2>u ) | L« = Xl , = x 2 '_ 

F (t + u 1; t + u 2 , (l« ^ - L$) + x a , {L^ U2>U - Lg) + x 2 ) ] . (5.8) 
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Since the LRB components are L-independent, and since the L-laws of the LRB components 
are known, the probability densities of the increments in the above equation are known. 
We thus have: 



L 



L 



(2) 
J t+u 2 ,U 



(2) 



exp 



.11 
2ui 



My 2 >o} mU2 _i 

— — y 2 exp(-y 2 )dy 2 , 

r[mu 2 ] 



(5.9) 



(5.10) 



where m > and T[x] is the gamma function. In order to work out an explicit example, we 
need to specify F(t,yi,y 2 ) and the weight function w(t, u 1 ,u 2 ). We choose the following 
example: 



F (t + u 1 ,t + U2,yi + x 1 ,y 2 + x 2 ) = exp [a (y 1 + x l ) -c (y 2 + x 2 )] , 
w(t, u 1; u 2 ) = exp {~\a\t + u x )) (c + 



(5.11) 
(5.12) 



where a e [—00,00), c > are constants. Then we insert (5.9) and (5.10) together with 
( |5.1ip and ( |5.12| ) in ( f577] ) and calculate the integrals over Ui and iz 2 . The result is: 

7 t L r = (17 - T) 2 (c + l) mt exp (a - c L ( t J - | a 2 t) . (5.13) 

The two-factor pricing kernel, jointly driven by a Brownian random bridge and a gamma 
random bridge, is thus given by a formula similar to (2.20) where {£ t } is the reciprocal of 
(5.2) while Y^ t , Y" L t , and Y" L are deduced from d5.13D. For this class of models, the bond 
price process can be written as follows: 



Pa 



P 0T + b{T)A\ 
P 0t + Kt)A\ 



where 



K0 = 



([/-t) 2 /i(0 



i+/i(0)t/ 2 ' 

A\ = {c + IT 1 exp [a L$ - c - \ a 2 t) - 1, 



(5.14) 

(5.15) 
(5.16) 



for < t < T < U. It is straightforward to show that {A^} is an ({^J,L)-martingale for 
t e [0, [/). Then the process {£ t A^} is an ({^" t },P)-martingale, and the bond price process 
(5.14) has a representation in terms of ({J^}, P)-martingales, that is 



P tT = 



0T*-t 



P 0t * t + b(0A*' 



(5.17) 



where A F t = £ t A^. Such models might be regarded as belonging to the finite-time equiva- 
lence class of Flesaker & Hughston (1996) bond price models. We note that pricing kernel 
models over infinite time, as in Flesaker & Hughston (1996), must be potentials of class D, 
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see Rogers (1997), Meyer (1962). In finite time, pricing kernel processes merely need to 
be positive ({J r t },P)-supermartingales to ensure non-negative interest rates. Furthermore, 
the martingale processes underlying the price models arise endogenously from the pricing 



kernel structure ( |2.1| ). Formula ( |2.1| ) can be viewed as a "machine" that implicitly produces 
martingales for pricing formulae with a rational form. Bermin (2012) revisits the Flesaker- 
Hughston approach to bond pricing and shows how yield curves may be inverted for any 
short rate process consistent with bond price processes that have an exponentially-affine 
structure. 

In Section [3] explicit pricing models are derived, and one may ask at this point what 
the connection is between those pricing models and the ones specified in (5.141. It turns 
out, as expected, that the link is a change of probability measure. Let us consider an 
({J ? t },L)-martingale {A^} and an ({J^J,M)-martingale {A^f}. Furthermore, we introduce 
an ({^"f}, L)-density-martingale {^ t }o<t<u tnat changes the probability measure L to the 
equivalent measure M. Then we set A 1 ^ = r] t A^f, and finally observe, for < s < t < U, 
that 



E L [A^ I JF ] = E L [q t Af | Jf] = r? s E M [a* | Jf ] = r? s Af 
This type of relation is also what connects ( [5.14 ) and ( |5.17[ ) . 



(5.18) 



A useful formula. The fact that an LRB has the law of its generating Levy process 
under the "Levy probability measure" L can be exploited to derive a rather useful formula 
involving the characteristic function of a Levy process and its Fourier transform. This for- 
mula might be particularly useful for numerical implementations of option pricing formulae 
based on Fast Fourier Transform techniques, see, e.g., Carr & Madan (1999). We consider 
again the L-conditional expectation 



E L [F(t + u,L t+u>u )\L tU ] , 
where {L tU } is an LRB. We specify F(t, x) by its Fourier transform F(t,y), that is 



(5.19) 



F(t,x) 



exp(-ixy)F(t,y)dy, 



(5.20) 



2 



where F(t,y) is selected such that F(t,x) is positive and integrable. Then we have, 



E L [F(t + u,L t+u>u )\L tU ] =E L 



exp(-iy L t+ uU ) F(t + u, y) dy 



nv 



(5.21) 



Assuming that Fubini's Theorem is herewith satisfied, we swap the expectation with the 
integral, and obtain 

- 

E L [F(t + u, L t+U}U ) | L tu ] = E L [exp (-iyi t+U;[7 ) | L tu ] F(t + u, y) dy. (5.22) 



16 



Since the LRB has the L-law of the underlying Levy process for t e [0, U), we can cal- 
culate the conditional expectation by recalling that the increments of a Levy process are 
independent and stationary: 



E L [exp (-iyL t+u>0 ) \ L tU ] = exp (-iyL tU ) E L [exp(-iyL uU )] 



(5.23) 



The expectation on the right-hand-side of the equation above is the generating function of 
a Levy process. We denote the characteristic function of a Levy process by \Ky), and thus 
write 

E L [exp (-iyL uU ) ] = exp [-u *(y )] , (5.24) 
for u e [0, U). This leads to 



E L [F(t + u,L t+u>u )\L tU 

and hence to the useful formula 

r u-t 



exp [-iyL tU - u *(u)] F(t + u, y) dy, (5.25) 



■ tT 



W 



T-t 
U-t 



T-t 



(T,u-r + t)E L [F(t + u,L t+UfU )\X tu ] du, 
w(T,u-T + t) exp [-iyL tU -u*(u)]F(t + u,y)dydu. 



(5.26) 



Expression (5.26) is valid also in the multi-factor case, in which the LRB {L tU } is a multi- 
dimensional vector. The elements of the LRB vector may be dependent through their termi- 
nal marginal laws as considered at the beginning of this section. The model (5.13 1 may be 
derived as a special case of the formula (5.26). 

We conclude this section by producing multi-dimensional and multi-factor asset price 
processes, and note that for the rest of this section the process {yQ} that follows is defined 
under the L-measure. For i = 1, 2, . . . , n, let 



tT ~ 



r U-t 



T-t 



E L [F f (t + u, L t+UtU ) | L tu ] Wi (T,u-T + t) du, 



(5.27) 



where Fj(t,x) is a positive and integrable function, and w ; (t,u) is a weight function. We 
emphasize that the Markov process {L tU } Q < t may be multi-dimensional. Then, the follow- 
ing is a multi-dimensional and multi-factor model for the bond price: 



P tT = 



where 



Por + sr = iy i (r)(yff-y c ?) 



w 

oo 



(5.28) 



(5.29) 
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These models can be extended further. Let us assume that the multivariate LRBs generat- 
ing the market filtration {3? t } satisfy Proposition 5.2 Then, the product of an QJ^^L)- 
supermartingale is again an ({^" t },L)-supermartingale. This leads us to the construction of 
higher-order asset pricing models, hereunder applied to the pricing of bonds. We generalise 
the bond price model ( |5.14| ): For < t < T and N e N, 



P tT = 



N 

i=i 7 Vr 

CO' 
t 



where 



A 



(0 

tT 



b« 



(T)Af ini 



In, 



1,-1 



Po t + Sf„A«' 
n;-0;-l) 

s 

Jn^— (m^— 1) Jn,-- 

,(0 



,0 



l(0 ^^(i.jiij-Cmj-l)) 



b w (r)A ( t i,j " i " l) 

Jn^ - 1 = jn i -2 + 1 Jn; = Jn i -1 + 1 



2 bf(T)4' ki) (5.31) 



(5.30) 



for n ; > m ; e N. By setting T = t in (5.31) one obtains {A tt }. The deterministic 
functions bf{t), bf{t), b®(t) are such that the product bf(t)b%\t) ■ ■ ■ bp it) is non- 
negative and non-increasing. The processes {A t } < t <u and {A^}o<t<u are ({^tl^) - 



martingales, and these are L-independent for ; ^ k. For instance, for N 
m 1 = 2, one obtains 



P tT — 



P 0T + b 133 mA™A™A™ 
P 0t + b lM A?A?A?' 



1. Hi = 3 and 



(5.32) 



where b 123 (t) = bi(r) b 2 (t) b 3 (t) for < t < T. For AT = 1, n x = 3 and m a = 1, we have 



P tT = 



P or + b 12 {T)A { ? A ( t 2) + b^m^A™ + b 23 (T)A ( t 2) A ( t 3) 
Pot + biaCOA? 3 ^ + bMA^A™ + b 23 (t)A ( t 2) A ( t 3) ' 



(5.33) 



where by(t) = bj(t)b k (t) for j ; ^ fc and < t < T. In order to lighten the notation, the 
i-index is suppressed in (5.32 1 and (5.33 1 since we have only one type of higher-order term 
in the sum over t. For N = 2, n 1 = 3 and m 1 = 2, n 2 = 3 and m 2 = 1, a combination of 
third-order and second-order models is obtained: 



P tT = 

P 0T + b^jm^A^A^ + b^jmf^A^ + b™(T)A?>VA™ + bg^m^Ap 

P 0t + bScOA^A?- 23 Ap» + b^WAp^ + blfWAF^p 5 + bg(tM^A? 3) 

(5.34) 

We note that the construction of higher-order pricing formulae is not limited to models 
driven by LRBs. The pricing kernel model ( |2.1[ ) can be used to construct higher-order price 
models driven by other Markov processes. Higher-order models gain in importance when 
considering general asset pricing including dependences across several types of assets. 
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6 General asset pricing in finite time 



The pricing formula ( TT7TT ) states that the no-arbitrage price process of an asset is an {{& t }, P)- 
martingale. There are several ways of constructing ({^ t },P)-martingales — we consider 
however a natural method within the framework developed thus far. For some fixed T , we 
denote by {m tT } < t < r an ({J^}, P)-martingale, and write 



n t S tT = m tT . 

Definition 6.1. Let {Z tr } < t < r<l/ be defined by 

c u-t 



Z tT = 



E F [G (t + u,X t+u ) \X t ] V(r, u-T + t)du, 



(6.1) 



(6.2) 



T-t 



where the deterministic function G(t,x) is measurable, and \p(t,u) is a deterministic function 
with the property *p(t,u — s) = \p(t —s,u)fors < t Au. 

Proposition 6.1. Let < t < T < U, and let g a (T) and gi(T) be deterministic functions. 
Then, for each fixed T, 

m tT = g Q {T) + gl {T)Z tT (6.3) 

is an ,F)-martingale. 

Proof. This proposition can be proven by following the steps in the proof of Proposition 
2.2 in Akahori & Macrina (2012). We observe that \p(t,u-s) = xp(t-s,u) implies ip(t, u) = 
^(t + u). □ 

We now have the necessary ingredients in order to propose the following class of asset 
price models. 



Proposition 6.2. Let {n t } and {m tT } be the processes (2.12) and (6.3), respectively. Then 
the asset price model (6.1) takes the form 

Sqt +z(T) (Z tT -Z ot ) 



S tT — 



Pot + y(t)(Y tt -Y 0t ) ' 



(6.4) 



where z{T) = gi(T)/7r . 



Proof. The relation (6.41 follows from (6.1) by inserting (2.12) and (6.3). The degree 



of freedom go(T) can be calibrated to the asset price S 0T at time via 

go(T) = S 0T n - g 1 (r)Z or , 



(6.5) 



where n - 1 + /i(0)Yo . □ 
As at the beginning of Section [3j we next derive general asset price models for which 
explicit expressions for derivatives can be computed. In doing so, we implicitly apply the 



results in Proposition 2.1 
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Lemma 6.1. Let {M t } < t<L/ be the {& t }-adapted density martingale inducing a change-of- 
measure from P to M. Let {A t }q<V<u ^ e ({^"t},M)-marttngaZe5, and consider a pricing kernel 
model of the form (3.1 K Furthermore, let 

Sj-j- 



s 0T + bxCrjAg 3 

P or + &2(nA ( r 2) ' 



(6.6) 



where b;(T), i — 1,2, are deterministic functions, and b 2 (t) is non-negative and non-increasing. 
Then the associated asset price process {S tT } is given by 



S 0T + b x (T)4 



(i) 



P 0t + b 2 (t)A { P ' 

Proof. The expression (6.7) is obtained by computing the expectation in 

1 



71, 



■E r [7i r S rr |«^ t ] 



(6.7) 



(6.8) 



where, by use of {M t }, the measure P is changed to M in order to exploit the martingale 



property of {A t } 



(0li=l,2 



t >0<t<U 



under ' 



□ 



Dynamical equation of the asset price process. We consider the price process given 
in (6.7) and assume that it is adapt ed to a market filtration generated by two Brownian 
random bridges {L^J}, i = 1,2, c.f. (3.9). We then follow the calculations in Section^to 
deduce the subsequent result: 

Proposition 6.3. Let {^ t } be jointly generated by {L^}, i = 1,2. Let the price process of an 



asset be of the form (6.7) where {A^o^^y satisfies 

dA ( t ° = v« (dW t F(i) + #«dt) 



for i = l,2. The process {v^} is \}-adapted, and 



C0_ u i 



U 



■ E 



dW ( 



p « = dL co 



hi" 



uw 



X 



(0 



•t —tU tj 

where u i is constant. The dynamical equation of such a price process is given by 

dS tT 



(6.9) 



dt 



where 



r t = 



P 0t + b 2 (t)A\ 



(2) 



P 0t + b 2 (t)A? y 



= (r t + A t S tr )dt + S tr dW t r , 



(6.10) 



Pij v t 



(2) 



P 0t +b 2 (t)A\- 



12! 



V 



(2) _ v (2) 



b 2 (0 
P 0t +b 2 (t)4 2) J 



12) 



V P 0t +b 2 (t)A ( t 2 V 



The process Wf = (W t P ,W t P ( ' 2 ' 1 ) is a two-dimensional ({& t },F)-Brownian motion where 
dW? (0 dW, 



p 0') _ 



p ; ■ dt/or i # and dW t P (0 dW t r UJ = dt for i = j. 



r* W _ . 



20 



As in Lemma 



4.1 



{V( }, i = 1,2, are determined by the specific choice of the models 



at the base of the processes {A t }, i = 1,2, and thus are model-specific. Furthermore, the 



dynamics (6.10) can be represented in the form so as to obtain the risk-neutral dynamical 
equation of the asset price process {S tT } < t < T . We have: 



dS tl 



= r t dt + £ tr dW7, 



(6.11) 



'tT 



where dW f 



dwf + A f dt is the risk-neutral Brownian motion defined in terms of the 



't — ""t 1 "t 1 

P-Brownian motion {W t F } and the market price of risk process {AJ. The solution to the 
stochastic differential equation (6.11 ) has the familiar Q-log-normal form 



$tT — Sot ex P 



/* 2 S sT 



)ds + 



r t 



X tT dW* 



(6.12) 



How one may obtain "finite-time Black-Scholes-type models" from ( 6.12 ) can be deduced 
by consulting Brody et al. (2008), Section 9. In addition, appropriate choices for the func- 
tions F(x), G(x), and the related weight functions will need to be made. 



Example. Let the market filtration be generat ed by two Brownian random bridges 



i = 1,2. Let us consider the quadratic model (3.11 1 for the pricing kernel modelled 



(2) 

by \A\ } driven by {L tV '}, and the exponential quadratic model (3.15 ) for the equity com- 



ponent modelled by {A^p} and driven by {L^j }. In such a case, one obtains an asset price 
process of the form 



'or 



+ 



Sf7 



4l/[l+(l/12) ACOjt/ 3 ] 



V^l - t/U exp 



r W 2 



2(1/- 



Pot + 



w[i+an2)fi(P)u 3 ] 



U t (2) 2 
(U-t) 2 tu 



t 

U-t 



In the case that {L tu } is multidimensional, that is, L tU = {L^ , L?}, L^}, the asset with 
price process ( |6.13| ) is traded in an incomplete market. 



(6.13) 



■Mi 



Example. Let the market filtration be generated by four LRBs: a stable- 1/2 random 
bridge {L^ }, two gamma random bridges {L^ } and {L^J }, and a Brownian random bridge 

{Ljy }• For details about 1/2-stable random bridges, we refer to Hoyle et al. (2010). The 
four random bridges are assumed to be independent under the "Levy probability measure" 



L. Next we make use of Definition |6.1| and Proposition 6.2 It is convenient to take the 



expectations under an auxiliary probability measure as presented in Proposition 2.1 Given 
that in this example the market filtration is generated by four LRBs, we choose to compute 
the expectations under the "Levy probability measure" L, under which the LRBs have the 
law of the generating underlying Levy process. For the function G(t + u,X t+u ) in Definition 



(6.1 1, we set 



G [t + u 1} t + u 2 , L% Ui>u , lfl Ui>u ) = exp [-KL^ + cL^) (6.14) 
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where k > 0, < c < 1. The Laplace transform of the stable-1/2 subordinator {L t }, which 
is a Levy process, is 

E [exp (— x:L t )] = exp I — — t I (6.15) 

where a > is the activity parameter that features in the subordinator's density function 

at 



p t (y) = My >o}- 



a 2 t 2 ' 



'2ny 



3/2 



exp 



2y 



(6.16) 



We now calculate (6.2) with the specification (6.14) and by setting 



ip(t + u) = ip(t, u l} u 2 ) = exp 
where m > 0. The result is: 



Z tT = ([/ - Tf exp ( -kL^ + ^ t ) (1 - c) mt exp (cl^) 



(t + u a ) (1-c) 



m(t+u 2 ) 
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(6.17) 



(6.18) 



In this case, the price process {S tT } given in (6.4) can be written in the form (6.7). For the 
denominator, we choose a discount factor of the kind (5.14). We obtain the following price 
process for, e.g., equity: 



S 0T + b x (T)4 



(i) 



P 0t + b 2 (t)Af ) ' 



(6.19) 



where, for r\ > 0, a e [-oo, oo), and q > 0, 



A { p = (1 - c) mt exp f -kL$ + ^ 



A { ? = (t, + 1)^ exp I -tjL^ j + aL^ - -a^t I - 1, b 2 (t) = 



(3) 



(4) 



([/-r) 2 gl (r) 
i + /i(0)t/ 2 ' 

l+ACOt/ 2 ■ 



(6.20) 



Another example could be constructed by generating equity models driven by "VG ran- 
dom bridges" and quoted in units of the natural numeraire at the basis of the model (5.14). 
In the case that the economic factors modelled by the random variables V^L are dependent, 
one obtains simple dependence structures between the dynamics of the equity and the as- 
sociated discount bond system that determines the discount rate in the financial market. 
The more advanced dependence models introduced in the next section can also be applied 
to model more complex interactions between different segments of a financial market. 



7 Spiralling debt and its impact on global bond markets 

We now address in more detail the pricing of sovereign bonds. Even though the majority of 
sovereign bonds pay coupons, we focus on discount bonds, for convenience. This simplifi- 
cation does not affect the view taken or the problem we intend to tackle here. The emphasis 
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is shifted on the value of a sovereign bond that should reflect the level of economic health 
of the issuing country. News regarding the bond market over the last few years has con- 
stantly reminded us that investors frequently balance the capability of a sovereign economy 
to grow vis-a-vis the amount of accumulated debt held at any one time. We choose this 
point of view, and wonder how to construct asset pricing models, which take into account 
at least some of this perspective. 

We consider a simple model for the economic structure of a country. We assume that 
a central government has a source of income, for instance taxes and the revenues of state- 
owned companies. On the other hand it also has expenditures in order, for instance, to 
finance armed forces, public education, a public health system, and other welfare. While the 
difference between income and expenditures fluctuates over the course of time, we assume 
that it is unlikely, at least in a well-run and periodically well-assessed economy, that this 
difference spikes for the better or for the worse. If it were the case, then we might see an 
economy's growth rate move from 1% to 10% within a few months, or a drop in the growth 
rate by a similar amount in the same time span. It is more likely though that a central 
government has to step-in to cover huge unexpected losses due to, e.g., the unfolding of 
an international financial crisis, domestic or international wars, natural catastrophes, and 
other calamities hard to predict and with disastrous impact on the economic health of 
a country. So, in addition to the "structural" income and outcomes of an economy, we 
consider the accumulation of significant debt due to severe losses, which may very well 
make the level of financial stress of an economy "jump". We model the structural part of 
the various cash flows of an economy by a Brownian random bridge 

L&WtXtf + fcu, (7.1) 

where X v may represent the economic wealth of a country at a future time U. We model 
the spiralling cumulative debt amassed by a sovereign country in the time interval [0, U] 
by a gamma random bridge {L tl j}. The random total (extraordinary) debt amassed by 

(2) (2) 

time U is modelled by Ly(j =X U . Its distribution can be arbitrarily specified. We imagine 
that the "structural" or "non-crisis" balance L^y is dependent on the total debt (losses) 

accumulated by time U. For instance, a sovereign government may decide at time U to 
make substantial cuts to the expenditures for public welfare if the amount of "extraordinary 
losses (debt)" will have spiralled by time U beyond what is perceived to be manageable. 
Therefore Lyjl and are assumed to have a joint marginal distribution, and this means 
that we are in the same modelling environment as in Section |5j 

The bond pricing model presented next is one of the simplest, though still rich enough 
to capture the desiderata within this discussion. One can of course choose to develop more 



sophisticated models. We choose a class of bond price models similar to (5.141, and follow 



the steps from (5.71 to (5.16) with one minor change in equation ( 5.11] ). We consider 



F{t + u x ,t + u 2 ,y l + x l ,y 2 + x 2 ) = exp (-a(y x + x x ) + c{y 2 + x 2 )) , 



a 2 



w (t,u 1 ,u 2 ) = exp[--(t + u 1 )Ul-cT (t+U2) . (7.2) 
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The reason for the change in the sign is that this way the losses will be recognised as 
downward jumps in the time series of the bond price. We emphasize that expectations 
are computed under the L-measure, under which LRBs with joint terminal marginals are 
nevertheless independent and inherit the law of the generating Levy processes. The bond 
price is then given by 

_P QT + b(T)A\ 

PtT ~ P 0t + Kt)A\> (7 - 3) 

where, for < t < T < U, a > 0, < c < 1, m > 0, we have 

W' lJ . A- = a-crexp(-aI«)-^ t + clS)-l. (7.4) 

Dependence in international markets. The effects of losses getting out of control 
are not confined to ones domestic economy. Especially in a global financial market, the 
deterioration of an economy's health exposes, for instance, foreign creditors holding debt 
of the distressed economy to higher credit risk. Bond markets are global "debt networks" 
linking several national economies with one another. The result of such network might be 
"contagion": an ailing economy may severely damage creditors which, through financial 
exposure, can be affected by spiralling losses experienced by the debtor. For instance, a 
foreign investor may see their investments in foreign bonds significantly devalued if the 
bond price declines due to unexpected losses, or out-of-control debt management. The 
foreign investor's loss may be commensurate with the percentage investment in an ailing 
economy compared with their total financial exposure and reserves gained through income. 
Here size matters, of course, and the discussions about the magnitude of the Greek versus 
the Italian debt impacting on the Eurozone or world economy come to ones mind. The 
next example aims at illustrating how contagion effects can be modelled in the present 
asset pricing framework. We introduce a linear combination of cumulative random bridge 
processes defined by 

^E^lS, (7.5) 

i 

where {L^} are, e.g., gamma random bridges with joint terminal distribution and gen- 
erated by independent gamma processes. The weight parameter measures the level 
of exposure to each cumulative process {if^}- Since the linear combination may not be 

the same for any creditor exposed to the pool further freedom is given through j- 

indexing the exposed entity (sovereign state, private company, etc.). For instance Country 
A may have a financial exposure of 15% to Country X and 8% to Country Y. The percentage 
exposures, which can be collected form various financial intelligence organizations, can be 
used to determine the weights for Country A. On the other hand, Country B maybe exposed 
by 34% to Country Ys economic performance and by 65% to Country 2Ts. One sees that, 
even though Country A and B may not hold any of each others financial assets, they are 
linked to each other through the common exposure to Counrty Ys economy, albeit to dif- 
ferent levels. The bonds of Country A and Country B are expected to both show the impact 
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of their respective exposures to Country Y. The bond price processes of Country A then may 
have the form fl7.3[ ) where 

([/-t)" +1 /i(0 



b(t) = 



1+A(0)I7 



^ * = n t 1 - w i j) ) h ^ z s) - 1 > 



1=1 



(7.6) 



and < < 1, m; > 0. The level of financial exposure to a specific economy can be mea- 
sured, at least in part, in terms of debt instruments held. Hereafter, we look at a simulation 
of contagion due to exposure to sovereign debt of two foreign countries. In particular, we 
suppose that Germany and France are exposed to the Spanish and Italian economic envi- 
ronment. The level of exposures w- , governing in part the levels of dependence among 
the debt holders, are selected as follows: 

GER FRA ESP ITA 
GER 10 
FRA 10 
ESP 0.57 0.47 1 
ITA 0.49 0.25 1 

In this basic illustration, Germany has no exposure to the French economy, however it has 
exposure levels 0.57 and 0.49 to Spain and Italy, respectively. The numbers listed in the 
table are not normalised. Dependence among the four economies is also subject to the joint 
marginal distribution of the LRBs underlying the dynamics of the yield processes. The closer 
the time horizon U, the more the joint marginal distribution of the multivariate random 
variable L UL i will govern the dynamics of the dependent yield processes. The simulations 
that follow are based on bond price models of the form (7.3) where £>(t) and {A^r} are 



specified in ( 7.6 1 
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Figure 1 : Simulation of the yield process of the two-year-maturity bonds issued by Germany, 
France, Italy, and Spain. 



In Figure [TJ we see the impact of spiralling debt on the yield process of sovereign 
bonds. The significant jumps in the yield processes of Spain and Italy (first and second 
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trajectories from above) are due to unexpected losses (e.g. bank bail-outs, natural disasters) 
over relatively short periods of time. These spikes then have repercussions on the yield 
processes of France and Germany to various degrees of severity. The level of repercussion 
on each exposed country is relative to the health of the own economy. In the simulation 
above, we see that France's yield process (third trajectory from above) needs more time 
than Germany's yield trajectory to recover from the shocks. The implication is that although 
Germany has a higher exposure level to Spain's and Italy's finances, it has a more robust 
domestic economy — with, e.g., higher growth rate — than France. Thus, Germany is in a 
position to better weather foreign economic shocks. Contagion effects, due to increased 
economic stress, are also observed in the behaviour of the spread process when comparing 
the performance of bonds issued by different sovereign states. 




2010 2011 2012 



Figure 2: Simulation of the price spread process for the two-year-maturity bonds issued by 
France, Italy, and Spain compared with the two-year bond issued by Germany. 

In Figure [2] we plot the difference between the yield processes of France, Spain, and 
Italy when compared with the yield process of the German sovereign bond. We observe the 
widening of the spread level for Spain and Italy, due to the two upward jumps in economic 
deficit. France keeps the evolution of the spread between the yield of its bond and the one 
by Germany in check even though it takes the hit from the exposures to the Spanish and 
Italian economies. 
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